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We present a review of the Quantum Monte Carlo results on the magnetic, charge and single-particle 
excitations as well as the pairing correlations of the two-dimensional Hubbard model. We are particularly 
interested in how these quantities are related to each other in the metallic state that forms near half- 
filling. These results are useful in gaining a better understanding of the low energy excitations of the 
superconducting cuprates. 



^ ■ 1 Introduction 

—5 ; 

. Superconducting cuprates have many anomalous physical properties Q . The parent compounds 

are antiferromagnetic (AF) charge-transfer insulators. Upon finite doping long-range order is 
destroyed, and a new strongly-correlated metallic state with unusual transport properties forms. 
In this state there are strong short-range AF correlations. Angular resolved photoemission ex- 
periments show that the single-particle properties are heavily damped and strongly renormalized. 
Charge dynamics probed by various transport measurements are also unusual. Most importantly, 
£N1 . a large number of experiments point out that the order parameter of the superconducting state 

' has d x 2_ y 2 symmetry [EL 

Perhaps, the simplest model that has similar electronic properties is the two-dimensional 
single-band Hubbard model given by 



> 



o 



H=-t^2 (4s c js + c] s Cis) + U ^n lT n 4 . (1) 
g ■ (ij>,» « 

Here c is creates an electron of spin s on site i, in the first term the sum is over near-neighbor sites, 
and n is — c| s Cj s is the occupation number for electrons with spin s on site i. The near-neighbor 
O ■ hopping matrix element is t, and the onsite Coulomb repulsion is U. 

Here, we will present Quantum Monte Carlo results on the magnetic, charge and single- 
particle excitations, and the pairing correlations of the Hubbard model in the intermediate cou- 
pling regime near half-filling for an 8 x 8 lattice. We will see that a strongly correlated metallic 
band develops upon doping the antiferromagnetic Mott-Hubbard insulator. This metallic state 
has anomalous electronic properties. For instance, the quasiparticle properties are strongly renor- 
malized. The system has strong short-range and low-frequency AF fluctuations, and the long 
wavelength charge response is enhanced. Furthermore, the system exhibits singlet d x 2_ y 2 pairing 
correlations There have been detailed Quantum Monte Carlo and exact diagonalization 

studies of these properties of the strongly correlated Hubbard or t—J models. Here, we are es- 
pecially interested in how these properties are related to each other. We will see that as the AF 
fluctuations grow, a narrow quasiparticle band, which has a large fermi surface, develops at the 
top of the lower Hubbard band. The AF fluctuations are also reflected in the effective particle- 
particle interaction, which is strongly repulsive at large momentum transfers. This feature of the 
effective interaction along with the large fermi surface causes d x 2_ y 2 pairing correlations 
We think that these numerical results are useful in gaining a better understanding of the nature 
of the low-energy electronic excitations in the superconducting cuprates. 
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Figure 1: Momentum dependence of (a) x(q, i<^m = 0) and (b) F(q) for U/t — 8 and (n) — 0.875. Here 
q is plotted along the path show in the inset of (b). 



2 Magnetic Fluctuations 

Using Quantum Monte Carlo simulations || we have measured the staggered magnetic suscep- 
tibility 



1 f 13 

a, J 



(2) 



Here mf(0) = c^c^ and m i+f (r) 



(0)e ffr , where m i+e (0) is the hermitian conjugate 



ofm+,(0). 

It is useful to have estimates of the characteristic length and energy scales of the AF corre- 
lations. For this purpose, in Fig. |l|(a) we first show x(Qj i^m — 0) versus q for (n) = 0.875 and 
U/t = 8. We see that as T is lowered below J ~ At 2 /U , strong AF correlations develop. Next, 
in Fig. 0(b) 

Imx(q,^) 

F(q) = hm (3) 

>o w 

is shown. Here, Im^fq, to) has been obtained by numerical analytic continuation JtJ of x(q, iw m ). 
It is well known that -F(q) is the quantity which determines the NMR T± response of the system. 
We see that F(q ~ (it, ir)) increases rapidly as T is lowered. The low frequency nature of the 
AF correlations are also seen in Fig. |^(a), where Im%(q = (ir, tt),w) versus ui is plotted. In the 
cuprates, the longitudinal relaxation rate Tj -1 of the planar Cu nuclei is dominated by the low 
frequency spin fluctuations near q = (tt, it). On the other hand, because of the form factor of the 
oxygen hyperfine coupling, the T^ 1 of the oxygen nuclei is mostly determined by the background 
of F(q) away from (ir,ir). 

In order to have an estimate of the characteristic length scales, in Fig. 2(b) we have plotted 
the magnitude of the magnetization correlation function |(m z (r)m z (0))| as a function of r along 
the (1,0) direction. We see that (m z (r)m z (0)) decays very rapidly, and the range of the AF 
correlations are of order a lattice spacing at these temperatures. 
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Figure 3: (a) Momentum dependence of x(«3j0)- (b) |(m z (r)m z (0))| versus r. These results are for 
(n) = 1.0 and U/t = 8. 
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Figure 4: Momentum dependence of n(q,0) for U/t = 8 at (a) (n) = 0.875 and (b) 1.0. 



These results demonstrate the existence of short-range and low-frequency AF correlations in 
the metallic state near half-filling. We note that while the AF correlations are short range, there 
is significant spectral weight under the peak at q = (tt, it) as seen in Fig. [I]. It is interesting to 
compare these AF fluctuations with those at half-filling, where the ground state has long-range 
AF order §. Fig. | shows results on %(q,0) and \{m z (r)m z (0))\ at half-filling. At T = 0.125t, 
the long-range order is clearly established on the 8x8 lattice. 



3 Charge Fluctuations 

In this section, we study the dynamics of the charge fluctuations. We will present results on the 
charge dynamics obtained from 

rf> 

n(q, iu m ) = \ dr e^ T (n(q, r)n(-q, 0)} , (4) 
Jo 

where n(q, 0) = ^ E P , CT c I+^ c p^- 

Figures ^(a) and (b) show Il(q, iui m — 0) versus q for (n) — 0.875 and 1.0. The dependence 
of II(q, 0) on U/t is shown in Fig. ||(a) for (n) = 0.875. Fig. ||(b) shows the temperature 
evolution of the charge fluctuation spectral weight ImII(q, u>) for q — (tt/4, 0) and (n) = 0.875. 
In these figures, we see that, on the average, the onsite Coulomb repulsion suppresses the charge 
fluctuations. At half-filling, II(q — > 0,0) vanishes as the Mott-Hubbard gap opens. However, for 
(n) = 0.875, the q — + part actually gets enhanced, as T is lowered. 

The q — > limit of II(q, 0) gives the compressibility n of the system. The low-temperature 
Monte Carlo studies |j| find that the enhancement of the long wavelength charge response near 
half-filling is due to the proximity of the system to a metal-insulator transition, and that, as the 
doping 8 is reduced, k diverges as <5 _1 . The enhancement of the long wavelength charge response 
has many experimental consequences. For instance, the long wavelength phonon modes of a real 
material can get softened by the enhanced II(q, w). 
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Figure 5: (a) Il(q, 0) versus q for various values of U/t at T/t = 0.5. (b) Imn(q,w) versus u at 
q = (7r/4, 0) and various temperatures for U/t = 8. These results are for (n) = 0.875. 



4 Single Particle Excitations 

In this section, we review numerical results on the single-particle excitations [|H|[n],[l2|] . We 
are especially interested in the effects of the spin and charge fluctuations on the single-particle 
excitations. The single-particle spectral weight and the density of states are given by 

A(p, uj) = — — Im G(p, iuj n — > uj + iS) (5) 

and 

iVH = l^A(p,c). (6) 
p 

In order to illustrate how the Coulomb repulsion affects the single-particle properties, in Fig. 
| we show A(p = (vr,0), uj) versus u for (n) = 0.875, T = 0.33t, and U/t = 4 and 8. Next, the 
temperature evolution of N(uj) is shown in Fig. |(a). The main features of the spectrum are 
lower and upper Hubbard bands, and a narrow metallic quasiparticle band at the top of the lower 
Hubbard band. Surprisingly, the quasiparticle band has strong T dependence. Comparing this 
figure with Figures ^ and ||(a) , we see that the temperature evolution of the quasiparticle band 
coincides with the development of the q ~ (ir,ir) magnetic and q ~ (0,0) charge fluctuations. 
Fig. |(b) shows N(lo) for U/t = 12 and T/t = 0.5, where the upper Hubbard band is further 
split from the lower Hubbard and quasiparticle bands. 

Results for N(u>) at half-filling are shown in Fig. pi At T — 0.5i, we observe a Mott-Hubbard 
pseudogap in the spectrum. At T — 0.125£, where long-range order is established on the 8x8 
lattice, a full single-particle gap is seen. At the moment, it is not completely understood how in 
the ground state of the sytem the metallic band, which is seen developing in Fig. 0(a), evolves to 
the N(ui) of the insulating state as the doping is reduced. At the temperatures that the Monte 
Carlo calculations are carried out, a metallic band forms at the top of the lower Hubbard band 
even for 3% doping. 
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Figure 9: Density of states at the fermi level N(0) versus the filling (n) for U/t = 8, and T/t = 0.5 
(solid circles) and f .0 (open circles). 
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Figure 10: t-matrix equation relating the reducible and the irreducible interactions in the BCS channel. 

Figure |^ shows the filling dependence of the density of states at the fermi level N(0) for 
U/t = 8. At T = 0.5t, the filling where the maximum of N(Q) occurs is near (n) — 0.85. 
However, the temperature evolution of the Monte Carlo data indicates that this critical filling 
will move closer to half- filling as T — > 0. In fact, for the simplest Hubbard model that we are 
using where only the near-neighbor hopping t and the onsite Coulomb repulsion U are taken into 
account, it might occur infinitesimally close to half-filling. Furthermore, simple additional terms 
in the Hubbard Hamiltonian, such as a next-near-neighbor hopping t', might also influence the 
position of the maximum of N(0). 

5 Pairing Correlations 

In this section, we will present results on the effective particle-particle interaction and the solution 
of the Bethe-Salpeter equation in the BCS channel. 

Using Quantum Monte Carlo simulations we have calculated the reducible particle-particle 
interaction T(p', —p',p, —p), which we will denote by T(p'\p) p3[ . Here, p stands for both p and 
iuj n . The reducible interaction is related to the irreducible interaction Ft through the i-matrix 
equation shown in Fig. |l0|. In the following, we will present results in the singlet channel where 
^s{p'\p) = |(r(p'|p) + F[—p'\p)). Figure |ll|(a) shows the momentum dependence of Ts(p'\p) 
and Tis(p'\p) at u> n = u> n > = irT, corresponding to an energy transfer uj m = 0, for U/t = 4 and 
(n) = 0.875. In this figure p is fixed at (w, 0) and q = p' p is swept along the (1, 1) direction. 
We observe that both Ts and Tig peak at large momentum transfers near q = (tt, 7r). The 
difference between Tg and Tis represents the effects of the repeated particle-particle scatterings 
in the BCS channel. Figure |Tll(b) shows rs(q, 0) versus q at various temperatures. We see that, 
as T is lowered, the q ~ (tt, tt) part of Fs gets enhanced, while the q ~ (0, 0) part decreases. 

Figure [l| compares the temperature evolution of Tisfq = p' P, i^m = 0) with that of 
x(q, 0) for U/t = 4. We find that rig(q, 0) shown in Fig. [L2|(a) can be fit exceptionally well with 
a phcnomcnological spin-fluctuation exchange form given by 

T l { q ) = U+\{ 9 Uf X ( q ), (7) 

where the renormalization coefficient g ~ 0.8. 

Tig can be used, along with the single-particle Green's function G(p,iui n ), to solve the 
Bethe-Salpeter equation 

KMp) = -^J2 t ^(p\p') g ^p') g i(~p')Mp')- (8) 

p' 
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Figure 11: (a) Momentum dependence of the particle-particle reducible and irreducible interactions in 
the singlet channel, Fs and Tis, at T/t = 0.25. (b) Temperature evolution of Fs(q, 0) versus q. These 
results are for U/t = 4 and (n) = 0.875. 
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The irreducible interaction Tis is attractive for the solutions which have positive eigenvalues. If 
the maximum eigenvalue reaches 1, then this signals a particle-particle instability of the system to 
a superconducting state. At the temperatures that the Monte Carlo calculations are carried out, 
the eigenvalues are much smaller than 1, and the system is far from a superconducting ground 
state, if one exists at all. Nevertheless, one can still study which pairing channels are favored 
as the AF fluctuations develop. We find that, in this regime, the singlet pairing correlations are 
strongest in the d x 2_ y 2 channel. The solid circles in Fig. [l3| show the d x 2_ y 2 eigenvalue as a 
function of temperature for U/t = 4 and (n) = 0.875. 

For U/t — 8, we have not been able to calculate the irreducible interaction because of nu- 
merical convergence problems, however we have calculated the reducible interaction and the 
Bethe-Salpeter eigenvalues. Figure |lj shows the temperature evolution of the reducible inter- 
action. Just as for U/t = 4, rg(q, iuj m = 0) at large momentum transfers is strongly repulsive, 
and it grows as T is lowered. Consequently, the leading singlet pairing correlations occur in the 
d x 2_ y 2 channel. The open circles in Fig. |l3| represent the d x 2_ y 2 eigenvalue for U/t = 8. 

Finally, we would like to note that the d x 2_ y 2 pairing correlations do not necessarily require a 
long AF correlation length, but simply large weight in T\s at momentum transfers near q = (n, tt). 
For instance, we have seen in Fig. |(b) that the AF correlation length for (n) = 0.875, U/t = 8 
and T/t — 0.33 is only of order one lattice spacing. 

6 Conclusions 

Numerical studies of the Hubbard model find that a correlated metallic band forms with unusual 
physical properties upon doping of the insulating half-filled system. Here, we have seen that 
in this metallic state the single-particle excitations are strongly modified by the many-body 
effects, the long-wavelength charge response is enhanced, and the system has strong short-range 
AF fluctuations. Also associated with this state are d x 2_ y 2 pairing correlations. 

These electronic properties are similar to those of the layered cuprates. Angular resolved 
photoemission experiments on the cuprates find that the quasiparticles are damped and renor- 
malized. NMR and inelastic neutron scattering experiments show the existence of short-range 
AF correlations in the superconducting samples. Furthermore, a large number of experiments 
indicate that the superconducting order parameter has d x 2_ y 2 symmetry. For these reasons, we 
think that simple strongly correlated models such as the two-dimensional Hubbard model are 
useful in understanding the unusual electronic properties of these materials. 
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